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I .  INTRODUCTION 


A  ballistic  algorithm  for  helicopter  fire  control  is  a  collection  of 
formulas  and  constants  which  in  conjunction  with  sensors  and  an  on-board 
computer  serve  to  predict  how  to  aim  the  aircraft  guns  and  rockets.  The 
ballistic  algorithm  for  the  2. 75  inch  rockets  as  developed  by  Dell  Helicopter 
*1  2 

Textron  (BHT)  *  is  based  on  a  modification  and  extension  of  equations 

3 

developed  by  Clark  of  BRL.  The  extension  by  BHT  consisted  primarily  of 
addition  of  equation  components  to  account  for  the  influence  of  downwash. 

The  BHT  modifications  included  alterations  made  in  the  structure  of  the  free 
flight  expression  for  "gravity  drop"  and  the  audition  of  ad  hoc  terms  to 
represent  the  effect  of  helicopter  velocity  normal  to  the  line  of  sight 
eferences  1  and  2  list  the  implemented  algorithm.  However,  there  seems  to  be 
nc  documentation  concerning  the  rationale  for  or  effect  on  accuracy  of  these 

4 

changes.  O'Bryon  of  BRL  performed  a  numerical  study  of  the  initial  modifica¬ 
tions  by  BHT,  Reference  1,  and  observed  significant  discrepancies  between 
sight  settings  as  predicted  by  BRL  versus  BHT  equations.  The  subsequent  BHT 
modifications  contained  in  Reference  2  differ  from  Reference  1  mostly  in  the 
readjustment  of  constants.  The  major  questions  regarding  the  BHT  equations 
and  modifications  concern  their  structure,  in  addition  to  the  specific  values 
of  certain  constants. 

One  of  the  major  changes  made  by  BHT  consisted  of  representing  grayity 
drop  by  an  expression  dependent  on  the  square  of  time  of  flight  (tA),  rather 

3 

than  a  power  series  in  range  as  derived  by  BRL.  A  conversation  between  the 
author  and  Mr.  Laird  Taylor,  formerly  of  BHT,  in  September  1980  indicates  that 
2 

the  t  term  was  empirically  based  but  was  thought  to  be  adequate. 

In  order  to  understand  the  effect  of  these  changes  and  determine  their 
validity  (or  not)  a  derivation  from  first  principles  is  made  herein.  This 

2 

derivation  will  show  how  a  t  term  can  be  justified  but  will  also  show  how  it 
yields  limited  accuracy  when  used  alone  (as  in  the  BHT  equations)  rather  than 
as  part  of  the  more  complete  ballistic  approximation  suggested  by  ballistic 
theory.  Simultaneously  the  development  will  shed  light  on  the  additional  ad 
hoc  terms  for  representing  the  influence  of  vertical  aircraft  velocity.  The 
development  will  also  show  the  omission  of  the  significant  influence  of  free 
flight  wind  effects  in  the  BHT  algorithm. 


*  . 

1.  Bell  Helicopter  Textron  Report  Ho.  209-099-520  (07/20/78),  " Modernized 
COBBA  Fire  Control  Computer  Mechanization  Equations,"  Rev.  B,  15  Aug  1978. 

2.  Ibid ,  Rev  E,  22  Aug  1979. 

3.  Clark,  David  L. ,  "A  Sight  Setting  Equation  for  Air  Launched  Ballistic 
Rockets,"  unpublished  BRL  manuscript. 

4.  O'Bryon,  James  F. ,  "Numerical  Comparison  of  BHT  Equation  # 209-099-520 
With  Comparable  BRL  Equation  of  the  MK  40  2.  75  Inch  Rocket,  Contained 
in  a  Itr  dtd  20  Sep  1578  to  PM  COBRA  (Mr.  Ryan),  SUBJECT:  Fire  Control 
Equation  Data  for  2.  75  Inch  Rockets  for  AH- IS  Modernized  COBRA . 
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II.  THE  IDEALIZED  FREE  FLIGHT  TRAJECTORY 


A  solution  for  an  idealized  free  flight  trajectory  will  be  developed 
parallel  to  that  obtained  by  the  author.  Ref.  5.  Here,  however,  the  coordin¬ 
ate  system  used  by  Clark  and  BHT  will  be  employed.  In  Figure  1  the  x  axis  is 
taken  to  be  along  the  initial  rocket  line.  The  z  axis  is  in  a  plane  formed  by 
the  rocket  line  and  the  gravity  vector,  z  is  normal  to  x  with  the  positive 
direction  downward. 


Figure  1.  Coordinate  system  used  in  developing  idealized  free  flight 
rocket  trajectory. 

By  proceeding  as  in  Ref.  5  the  differential  equations  governing  the  ideal¬ 
ized  trajectory  are  given  by 


x  ■  -n  (x  -  w  )  +  g 

A  A 


z  =  -  n  (x  -  w  )(z  -  wj  +  g 

A  L  l* 


(1) 
f  2) 


Eqs.  (1)  and  (2)  are  to  be  solved  subject  to  the  initial  conditions 


breaux,  Harold  J. ,  "A  Methodology  fop  the  Development  of  Fire  Control 
Equattons  for  Guns  and  Rockets  Fired  From  Aircraft,  •*'  April  1982  Draft 
BRL  Manuscript.  3  J 
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Here  the  subscript  b  denotes  burnout  conditions.  In  Eqs.  (1)  and  (2)  n  is  a 
lumped  parameter  given  by 

D  =  ird^  Cp  (M,  a)/8  m  ,  (4) 

d  is  the  rocket  diameter,  m  the  rocket  mass  and  Cp  (M,  a)  is  the  drag  coeffi¬ 
cient.  M  is  the  mach  number  and  a  is  the  angle  of  attack.  The  notation 
Cp(M, a)  is  used  to  indicate  that  Cp  is  evaluated  for  the  launch  mach  number  M 

and  launch  angle  of  attack  a.  Cp  is  thereby  assumed  constant  leading  to  n 

also  being  constant.  (NOTE:  The  complete  methodology  (not  discussed  here) 
provides  a  means  to  compensate  for  the  error  associated  with  this  assumption.) 
To  obtain  a  solution  to  Eqs.  (1)  and  (2)  one  can  proceed  as  in  Ref.  5.  Let 


C  X 

"  Xb 

-  w  (t 

X 

-  V 

(5) 

K  =  z 

-  Zb  • 

WZ  Ct 

-  V 

• 

(6) 

After  these 

transformations  Eqs. 

(1)  and  (2)  transform  to 

•  • 

C  “  - 

n  l2  + 

gx 

(7) 

•  • 

E  3  - 

•  • 

n  C  5 

-  gz  • 

(8) 

The  initial 

conditions  transform 

to 

£<tb) 

“  Vxb 

-  w  = 

X 

Vxb 

(9) 

‘<V 

=  0 

(10) 

Icy 

=  Vzb 

-  w  = 
z 

Vzb 

(ID 

scy 

=  0  , 

(12) 

To  obtain  the  solution  to  Eq.  (7)  first  introduce  the  transformation 


C 


[u(c) 3 1/2 


(13) 
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Inserting  Eq.  (13)  into  (7)  yields  the  transformed  differential  equation 


du/d;  =  -2n  u  +  2gx> 


(14) 


which  has  the  solution 


u  "  g  /n  +  u  2  e”2ri? 
6x'  1  xb 


(15) 


uxb  is  given  by 

*  o  2 

u  ,=  v  ,  -  g  /n  . 
xb  xb  sx ' 


(16) 


Since 


V  ,  >>  g  / n 

xb  sx'  1 


(17) 


hereafter  the  solution  wil.1  be  written 


u  =  g  /n  +  v2  e"2r|C 
6x'  xb 


(18) 


By  employing  Eq.  (18)  in  (13)  one  obtains  the  solution  for  t 

'-‘b  -/'  iv  *  vxb  «'2n'i‘1/2<i<  ■  nVr1/2  dp 

0  •'l 


where 


1  ■  VKb’  • 


(20) 


The  integral  in  Eq.  (19)  is  given  by 


t  -  t.  =  (nv  .  A172)"1  in  [pX1/2  +  (X  p2+l)1/2] 


b  "  '•"’xb 
Now  observe  that 


(21) 


In  [q  +  (q2+l) 1/2]  =  q-q3/5  +  3q5 


/ 40  -  .  . 


(19) 


(22) 


By  use  of  Eq.  (22)  in  (21)  one  obtains 


T 

t  -  tb  =  (nvxb)_1{(enC  -  1)  -  (X/6)(e3n?  -  1)}  . 

(23) 

.  4. 

In  seeking  a  solution  for  Eq.  (8)  we  let 

=  d5/dc  . 

(24) 

B  " 

One  has 

'* ' 

►* 

4pnr* 

II 

4P1* 

(25) 

i 

••  *  *  •  2  1  »2  * 

5  »  5  C  -  nS  C 

(26) 

leading  to  the  transformed  equation  for  (8) 

*1.  ** 

I  V 

C"  =  gz/C2  «  gz  e2nC  . 

(27) 

’ 

i' 

This  equation  can  be  integrated  directly  to  yield 

r. 

=  gz(e2nc  -  D/[2nv2b)+  c'(tb) 

(28) 

. 

►! 

where 

a 

C'<V  =  «/5lt  =  tb  ■  v*b/vxb 

(29) 

■4 

•  and 

V 

i  ' 

b‘* 

f- 

v  ,  =  V  .  -  w 
zb  zb  z 

(30) 

Lv 

h 

i*’ 

A  second  integral  yields 

?  =  gI(e2nVl)/[4n2v2b,  .  ^  /[2tw2b]  .  vJbt/vxb  . 

(31) 

ji 

If 

NOTE:  Here  gx  has  been  dropped. 

►  i 


i 
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This  equation  is  the  solution  for  the  normalized  drop  £  in  terms  of  normalized 
range  £.  To  change  the  dependence  of  £  from  c  to  t  one  can  proceed  as  follows 
Observe  that  in  Eq.  (23)  the  dominant  part  of  the  solution  can  be  written 


» -  -  owxbr  Vc-u  (32) 

since  X  is  small  for  launch  conditions  having  small  angles  with  respect  to 
horizontal.  This  can  be  viewed  as  a  first  order  solution  to  Eq.  (23).  To  in¬ 
corporate  the  X  dependence  let 

en?  =  1  +  y  (33) 

where 

r  -  nvxb(t-tb)  .  (34) 

Also  by  squaring  and  cubing  Eq.  (33) 

e2n5  -1  -  2y  +  Y2  (35) 

e3n5  -  1  =  3y  +  3y2  +  y3  .  (36) 

From  Eq.  (23)  one  has 

en?-l  .  y  +  (A/6)(e3n?-l)  =  y  +  (X/6) (3y  +  3y2  +  y3) .  (37) 

By  squaring  Eq.  (37)  one  obtains 

e  -1  =  2y  +  y  +  X(y  +  2y  )  +  term  of  order  (Xy  )  +  term  of  order  (X^). 

By  use  of  Eq.  (38)  (without  the  higher  order  terms  in  Eq.  (31)  one  obtains 

S  -  gz(2Yn2+X{y+2Y2})/[4n2Vx2]-gzi;/[2nv2b]+vzbC/vxb  .  (39) 

By  restoring  all  variables  to  their  original  form  and  observing  that 

vzb'  xb  zb  xb  z  xb  * 
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(40) 


Eq.  (39)  takes  the  form 


*■  s  *".(*-  »b  -  * 


xb'  ■  2nvxb  (t  ■  lb  ■  s',Vxb) 


+  ■;  g  (t  -  t,  )  +  - - 

4  &z  K  b'  4riv, 


2  gzHt-tb)  gzX(t-tb) 

- 4 -  +  Vzb  c/vxb 


(41) 


xb 


This  equation  provides  a  precise  solution  form  for  analyzing  the  deficiencies, 
omissions  and  correctness/o*  not  of  the  BHT  free  flight  components  of  the 
COBRA  rocket  ballistic  algorithm. 


III.  COMPONENT  ANALYSIS  OF  THE  GENERAL  SOLUTION 


z  is  the  so-called  drop  of  the  rocket  in  a  vertical  plane  normal  to  the 
initial  rocket  line;  zb  is  the  drop  which  occurs  during  the  xocket  boost  phase. 

The  term  wz  (t  -  tb  -c/vxb)  is  the  effect  caused  by  a  wind  normal  to  the 
rocket  line  in  free  flight. 

2 

Except  for  the  coefficient  of  1/4,  the  term  1/4  g^ft  -  tb)  corresponds 

to  the  free  fall  of  an  object  in  vacuum.  The  term  [gz/2nvxb] (t  -  tb  -  C/vxb) 

represents  the  change  in  free  fall  due  to  drag.  It  is  important  to  note  that 
as  drag  becomes  very  small 

t  -  tb  c/vxb  .  (42) 

Also  as  drag  approaches  zero  the  quantity  n  also  approaches  zero  so  that  the 
term  as  a  whole  becomes  indeterminate.  It  can  be  shown  that 


lim 

n-Hi 


I*  -  *b  -  5/v*b1  ■  "T  <* 


(43) 


This  result  is  obtained  by  expanding  the  exponential  in  Eq.  (32)  in  a  Taylor 
series  and  inserting  the  result  in  Eq.  (43).  Accordingly  as  n  approaches  zero 
the  sum  of  the  two  terms  approach  1/2  gz  (t  -  tb;  as  expected. 

2 

The  two  terms  1/4  gz*(t  -  tb)  and  [gzX(t-tb)l/4nvxb  are  interactions  be¬ 
tween  drop  in  -,he  z  and  x  directions.  For  flat  trajectories  that  are  nearly 
horizontal  these  terms  are  relatively  small.  Through  regression  analysis  it 

2 

has  been  determined  that  1/4  gz-\(t-tb)  is  the  more  significant  term  and  the 

remainder  can  be  dropped  for  the  2,75  inch  rocket  at  ranges  limited  to  6  km  or 
less. 
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The  term  Vzb  C/vxb  is  closely  related  to  the  burnout  angle  of  the  velocity 
vector.  Note  that 

c  a  X  -  -  wx  (t  ••  tb)  (44) 


and 


vxb  ■  Vxb  -  “x  • 

For  zero  wind  along  the  x  direction,  w^  =  0  and 
Vzb  ?/Vxb  *  <W  (X  -  V  s  6b  (x  -  xb> 


(45) 


(46) 


where  0^  is  the  velocity  burnout  angle.  The  availability  of  this  component 

analysis  now  makes  possible  a  precise  critique  of  the  BHT  free  flight  compon¬ 
ents. 


IV.  THE  BHT  FREE  FLIGHT  COMPONENT 

The  BHT  algorithm  is  designed  to  provide  the  angular  setting  for  engaging 
a  moving  target  at  a  range  R  having  a  velocity  Vt&r  along  the  line  of  sight 

r 

and  V  normal  to  the  line  of  sight.  Warhead  effectiveness  criteria  dictate 
X3r 

E 

that  the  rocket  warhead  functions  at  an  offset  point  given  by  the  increments 
AZ  and  AR  from  the  target  location.  The  R'  range  is  defined  to  be  the  range 
between  the  point  of  firing  and  the  desired  point  of  warhead  event.  The  latter 
point  is  affected  by  target  motion  and  offset.  Hence 


R'  =  R  +  V  t  +  AR  (47) 

tarr 

BHT  defines  an  angle  e  as  the  total  angular  drop  from  effective  launch  line 
to  impact  (warhead  event).  This  angle  is  given  by 

e  =  sin-1  ZBHT/R'  (48) 


ZgHT  is  approximated  by  the  equation 


ZBHT  "  cos  Ge^Zb  +  AZ  +  C  (,-V  1  +  VR’  “  V-Vtar  t  +  VAZ  1 


(49) 
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where 


* 


V 

gw 


V  tan  6 

gu 


(50) 


The  components  of  V ^  are  defined  as  follows:  u,  v,  w  are  body  axes  components 

for  the  helicopter,  u  is  forward  along  the  armament  datum  line  (ADL),  v  is 

through  the  right  wing  and  w  points  down  (normal  to  both  u  and  v) .  The  sub¬ 
script  g  refers  to  velocity  with  respect  to  ground,  then  resolved  in  this  u, 

v,  w  system.  0  is  the  angle  between  the  ADL  and  horizontal.  V  and  V 

gw  gu 

are  helicopter  velocity  components  resolved  as  described  above. 

By  comparing  Eq.  (49)  with  Eq.  (41)  and  in  light  of  the  component  analysis 
the  following  major  differences  are  evident. 

A.  The  BHT  equation  has  no  drag  term.  In  an  early  version  of  Ref.  1,  C 
was  set  to  g/2.  Later  it  was  adjusted  to  a  value  of  4.1  (=  .84  g/2). 

B.  There  is  no  compensation  for  wind  normal  to  the  launch  line  in  free 

? light .  "(Absence  o£  wz  (t  -  t^  -  c/v^)  or  equivalent  term) . 

C.  The  VAZ  t_  term  seems  to  be  nad  hoc"  in  that  it.  was  added  without  phy¬ 
sical  justification  or  fitted  correlation.  The  conjecture  here  is 
that  BHT  reasoned  (faultily?)  as  follows:  [The  vertical  component  of 
helicopter  velocity  is  given  approximately  by 


VAZ(vert)  *  V  '  Vgu  ta"  6  •  (51) 

This  velocity  at  launch  is  effectively  imparted  to  the  burnout  velocity 
and  the  decrement  in  velocity  remains  throughout  free  flight,  leading 
to  a  total  d  mnent  in  vertical  displacement  given  by  V^t.]  (End 
of  conjectur  Validated  ballistic  theory  suggests  that  the  effect 
of  any  such  velocity  component  is  the  so-called  "weathervaning"  effect. 
The  imparted  velocity  at  burnout  is  of  sign  opposite  to  the  effective 
wind  and  its  magnitude  is  multiplied  several  fold.  More  correctly  it 
should  enter  through  the  quantity  V  contained  in  Eq.  (41).  That 

equation  shows  that  the  general  solution  depends  on  the  initial  z 
velocity  (the  burnout  component  V  ^)  only  through  the  term  vz^  • 

Hence  a  z  component  of  velocity  such  as  V._  at  launch  in  inertial 

space  leads  to  a  z  component  V  at  burnout  which  as  noted  above  is 

^  AZ  D 

quite  different  from  V.7  due  to  weathervaning.  The  free  flight  effect 
of  VAZat  launch  leads  to  addition  of  a  term  given  by 


★  it 

VAzb 


** 


5/vxb  =  VAzb  *  Cc/v 


xb 


t) 


it  it 

VAzb  t  ^^vacuunP 


(52) 
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which  arises  directly  from  Eq.  (41).  Hence,  in  addition  to  the  dif- 
ference  between  and  V^z,  the  expected  effect  differs  by  the 
above  noted  range  ratio. 


In  any  case  the  effect  of  initial  z  velocity  can  only  be  obtained 
by  a  self  consistent  fitting  of  data  which  includes  weathervaning  in¬ 
fluence  in  addition  to  any  translational  effects.  The  ad  hoc  addi¬ 
tion  to  a  formula  of  such  a  translational  term  after  fitting  can 
only  be  done  correctly  if  the  data  did  not  contain  the  translational 
effect.  The  term  in  question  ijs  believed  to  be  a  significant  source 
of  error  in  the  COBRA  rocket  ballistics. 


D. 


There  is  no  compensation  in  the  time  of  flight  equation  for  range 
winds .  As  a  result  thereTs  no  compensation  in  the  sigh 


_  >ight  settings  for 

the  effect  of  range  winds  since  the  free  flight  influence  would 
correctly  enter  through  the  t  terms  in  Eqs.  (41)  and  (49).  The  audit 
trail  of  this  historical  omission  reveals  the  following.  The  source 
of  the  BHT  time  of  flight  equation  is  the  BRL  document,  Ref.  6.  That 
equation  was  developed  under  the  assumption  of  no  wind.  It^  did  not 
include  a  wind  correction.  That  equation  does  contain  a  correction 


for  non- zero  helicopter  velocity  V  .  That~correction  takes  the  form 


At 


V2(r-V3  Vvb 


(53) 


where  p  is  air  density,  V 


is  helicopter  speed,  and 


velocity.  As  is  a  fitting  constant. 


is  burnout 
A  telephone  conversation  with 


Laird  Taylor,  formerly  of  BHT,  reveals  the  following.  Taylor  indi¬ 
cated  that  the  time  of  flight  equation  was  fitted  to  trajectory  data 
which  included  significant  influence  of  wind.  The  above  term  was 


modified  by  replacing 


V  with  V  „  where 
a  au 


au 


V  -  w 
gu  u 


(54) 


BHT  apparently  believed  that  one  could  make  the  time  of  flight  equa¬ 
tion  applicable  to  wind  conditions  by  merely  replacing  the  ground 
speed  Va  with  the  air  speed  Vflu.  For  example,  BHT  in  Ref.  2,  page 

4-3  states  that  "For  rockets,  wind  is  implicitly  contained  within  the 
BRL  equation."  This  is  not  correct.  The  dominant  effect  of  wind  is 
corrected  for  by  the  general  solution  for  t  through  the  definition  of 
C  and  v  For  more  details  see,  Ref.  5.  The  omission  of  wind 

correction  for  time  of  flight  in  free  flight  is  a  major  source  of 
error  in  the  BHT  COBRA~~rocket~Fal  listic  algorithm.  (Note  for  the 
MK  66  with  submunition  warheads,  warhead  event  is  based  on  fuze  set¬ 
ting.  This  fuze  setting  can  be  in  error  by  about  .5  seconds  for  a  20 
knot  range  wind)  . 

letteA  ^ tom  S ft'L  to  USA  MZCOM  (Ma.  BeAgmn) ,  SUBJECT:  "Time  oi  F. tight 
Equation*  ion  2. 75  Inch  Rocket/'  22  Pec  1976. 


E.  The  BHT  equations  do  not  compensate  for  cross  winds  in  free  flight. 

This  observation,  follows  immediately  from  considerations  parallel 
to  the  above  discussion.  The  displacement  in  impact  prediction 
(error)  is  closely  approximated  by  w  (t-t.-  ?/v  .  )  where  w  is  the 
cross  wind.  y  x  y 

*  * 

F.  The  term  t  is  symmetrical  in  structure  to  V^z  t  and  is  also  a 

significant  source  of  error. 

G.  Through  conversations  with  BHT  personnel  it  is  apparent  that  the 
weathervaning  term  for  the  azimuth  equation  (for  the  MK  40)  was  field 
adjusted  to  "correct"  for  wind  effects  showing  a  bias.  In  the  ab¬ 
sence  of  a  predictor  for  the  free  flight  wind  effect  this  field  ad¬ 
justment  is  probably  in  error.  The  problem  is  compounded  due  to 
symmetry  since  the  adjustment  thereafter  (and  currently)  also  is  in¬ 
cluded  in  the  equation  for  the  elevation  angle. 


V.  NUMERICAL  ANALYSIS  OF  THE  DEFICIENCIES 

The  numerical  effect  of  the  deficiencies  related  to  drag  and  wind  effects 
in  free  flight  can  be  studied  by  correlation  of  trajectory  data.  Table  I  is 
a  listing  of  trajectory  conditions  for  an  extensive  file  of  MK  40  data  at  BRL. 
The  free  flight  portion  of  this  trajectory  data  can  be  correlated  with  Eq.  (41). 
This  equation  can  be  rewritten  in  the  form 


f  *  Vzb^Vxb"^V»h/,V»bHCz-zh^“W7^t“th''i::/V5fh^/5 


Al(vxb//Vxb^gz/,nvxb^t"tb  ~c/vxb*  +  A2(vxb/Vxb^gz^"  V 


(55) 


The  following  should  be  noted  in  this  equation.  The  quantity  vzt/vxtj  is  essen¬ 
tially  the  burnout  angle  of  the  velocity  vector.  The  quantity  vxb/Vxb  is 
nearly  unity  since 


vxb/Vxb  ■  <Vxb-\>/Vxb  “O 


V  ,  >>  w 

xb  x 


(56) 


The  quantity  (z-z^)/?  is  an  angle  associated  with  the  target  height  z.  The 

quantity  wz  (t-t^-c/v^)  is  a  drift  term  associated  with  vertical  wind.  The 

equation  is  dependent  on  range  winds  due  to  t  correctly  representing  the 
effect  on  time  of  flight  of  a  range  wind  and  on  the  definition  of  c,  i.e., 


c  =  X  -  Xb  -  wx  (t  -  tb) 


(57) 
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It  can  be  shown  that  the  BHT  rocket  algorithm  is  equivalent  to  the  above  equa¬ 
tion  with  all  wind  effects  absent  and  with  A^  =  0. 

To  study  the  effect  of  the  individual  deficiencies  the  data  referred  to 
above  was  correlated  under  various  assumptions.  In  Figure  1  the  correlation 
is  first  performed  with  all  wind  terms  present  but  only  the  t  term  for  the 
gravity  drop.  The  sliding  bias  is  evident.  The  single  term  cannot  do  the 
job  for  all  ranges.  Residuals  are  biased  positive  for  short  ranges  ind  nega¬ 
tive  for  long  ranges.  The  addition  of  the  second  term  in  Figure  2  r  moves  the 
bias  and  greatly  tightens  the  correlation.  Figure  3  shows  the  effeu;  of  the 
single  term  correlation  when  the  wind  components  of  the  model  are  removed  as 
in  the  case  with  the  BHT  algorithm.  The  maximum  error  increases  from  16  mils 
to  38  mils. 

2 

Figure  4  shows  the  effect  of  utilizing  the  coefficient  of  the  t  term  as 
is  currently  used  in  the  modernized  COBRA  for  the  MK  40.  The  correlation  (with¬ 
out  adjustment)  is  excellent  for  ranges  less  than  3  km  but  thereafter  a  growing 
bias  is  encountered.  At  6  km  the  bias  is  about  50  mils  with  a  maximum  error 
of  60  mils.  Figure  5  indicates  that  when  winds  are  also  included  the  bias  at 
6  km  is  about  60  mils  with  a  maximum  error  of  80  mils. 


The  deficiencies  cited  in  the  BHT  free  flight  model  can  be  overcome  by 
utilization  of  the  approximation  represented  by  Equation  (41)  restructured 
into  the  form  of  Equation  (55) .  The  correlation  obtained  with  the  three 
dominant  terms  is  shown  in  Figure  6,  The  RMS  and  maximum  error  are  reduced 
by  a  factor  of  17  from  that  depicted  in  Figure  5. 

It  should  be  noted  that  the  above  analysis  is  contingent  on  all  boost 
phase  effects  being  perfectly  represented  and  that  no  "synergistic"  interac¬ 
tions  arise  during  the  iterations.  As  shown  herein  there  are  additional  pro¬ 
blems  with  the  algorithm  due  to  poor  boost  phase  modeling.  This  of  course 
compounds  the  problem. 


VI.  AN  ANALYSIS  OF  PAST  ROCKET  FIRING  TESTS 


Conversations  with  test  directors  and  technical  personnel  associated 
with  past  firings  of  the  2.75  inch  rockets  from  the  COBRA,  AH-1S  indicate  that 
firings  show  no  obvious  bias.  This  is  in  contrast  to  the  inferences  made 
herein.  Some  shots  (means  of  a  group)  art  long  and  others  are  short  -  as  to 
be  expected  from  such  an  inaccurate  weapon  system.  Despite  these  past  cursory 
observations  the  analysis  performed  herein  and  the  audit  trail  of  modeling 
decisions  made  by  BHT  lend  insight  which  leads  to  the  discovery  of  a  large 
bias  in  ballistic  performance  heretofore  undetected.  As  will  be  shown  below 
this  bias  is  as  high  as  60  mils  in  elevation  at  a  range  of  6  kilometers. 

This  in  turn  leads  to  a  bias  in  rocket  impact  of  over  400  meters. 


Sections  III,  IV  and  V  showed  that  the  BHT  modeling  of  gravity  drop 
ignores  a  very  significant  effect  arising  from  drag.  Additionally,  a  constant 
in  the  weathervaning  term  was  empirically  adjusted  to  match  test  data.  Such 
an  adjustment  of  one  constant  could  only  produce  a  compromise  since  three 
significant  defects  exist.  The  third,  in  addition  to  the  drag  and  weather¬ 
vaning,  is  that  related  to  the  initial  condition  effect  for  velocity  components 
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normal  to  the  line  of  sight.  Intuition  suggests  that  the  empirical  adjusting 
of  the  weathervaning  constant  of  necessity  also  adjusted  for  the  mismodeling 
of  the  other  two  effects.  TTiis  implies  that  for  hover  firings,  wheve 
weathervaning  is  generally  insignificant,  test  firings  would  show  a  bias  due 
primarily  to  the  mismodeling  of  gravity  drop. 

Test  data  for  2.75  inch  MK  40  rocke^fir^.ngs  employing  the  BHT  algorithm 
are  contained  in  TECOM  reports  by  Andrese/  and  Sanborn8.  These  reports  convert 
miss  distance  to  an  equivalent  elevation  error.  As  a  result  of  insight  pro¬ 
vided  by  the  above  discussion  the  data  for  hover  firings  were  examined  in 
search  of  the  inferred  bias.  Test  data  are  superimposed  upon  the  elevation 
error  predictions  shown  in  Figure  S.  The  test  data  provide  conclusive  con¬ 
firmation  of  the  bias  inferred  from  the  theoretical  ballistic  analysis  employed 
in  the  critique.  The  bias  at  6  km  is  approximately  60  mils,  but  as  inferred 
by  the  analysis  the  actual  error  can  be  as  large  as  80  mils  due  to  the 
combination  of  effects  when  wind  is  also  included. 


VII.  CONCLUSIONS 

The  accuracy  of  the  Cobra  2.75  inch  rocket  systems  is  needlessly  reduced 
by  deficiencies  in  basic  equations  in  the  ballistic  algorithm.  A  question 
concerning  the  adequacy  of  changes  made  by  Bell  Helicopter  Textron  to  BRL 
provided  ballistic  equations  is  largely  answered  by  theoretical  analysis  and 
confirmation  utilizing  test  data.  Deficiencies  such  as  the  lack  of  a  drag 
component,  lack  of  wind  corrections  in  free  flight,  and  incorrect  compensation 
for  helicopter  velocity  normal  to  the  line  of  sight  can  lead  to  as  much  as 
80  mils  elevation  error  (£  600  meters  range  error)  at  maximum  range  fo-  the 
currently  implemented  MK  40.  Use  of  the  structure  of  the  existing  algorithm 
for  the  MK  66  would  needlessly  reduce  the  system  accuracy  for  that  rocket 
family .  A  means  for  correcting  most  of  the  deficiencies  is  available  and 
currently  being  implemented  through  the  analysis  efforts  recently  performed 
at  BRL  in  support  of  the  AAH  (Apache)  development. 


J.  Andrese,  Final  Independent  Evaluation  Report  for  the  Enhanced  Cobra 
Armament  Program  (ECAP)  of  the  AH-1S  Modernized  COBRA,  US  Army  Test  and 
Evaluation  Command,  APG,  MD,  Sep  1980. 

g 

J.  Sanborn,  Final  Report  of  First  Article  -  Preproduction  Test  (FA-PPT)  of 
the  AH- IS  Modernized  COBRA,  TECOM  Project  No.  4-AI-100-01S-01 7,  Oct  1980. 


TABLE  1.  MK  40  (M1S1)  TRAJECTORY  DATA  FILE 
(FILE  2755) 

DATA  BASE  OF  SIX  DEGREE  OF  FREEDOM  TRAJECTORY  CALCULATIONS 


VA 

(Knotc ) 

DOWNWASH  * 
VELOCITY 
(M/S) 

(MK  40  2 

P 

( X  of  Std) 

!.75  "  Trajectories) 

W1  w3 

(Knots)  (Knots) 

Range  Cross 

Wind  Wind 

DIVE 

Deg. 

ALT 

of 

Hel. 

M 

l.C. 

(Code) 

s  ,*45 

•  / 

0 

+102 

0 

10 

0 

15.2 

01 

-20 

0 

Std. 

10 

-10 

-5 

500 

02 

0 

13.716 

-102 

0 

0 

-10 

1000 

03 

15 

6. 058 

Std. 

0 

0 

0 

2000 

04 

30 

0 

+52 

-20 

-10 

+5 

15.2 

05 

50 

0 

-52 

20 

-20 

+10 

100 

06 

100 

0 

Std. 

-10 

0 

-10 

1000 

07 

150 

0 

Std. 

0 

20 

+10 

15.2 

08 

200 

0 

Std. 

0 

10 

0 

250 

09 

8 

10.058 

+  10% 

0 

-20 

10 

750 

10 

23 

3.200 

-10% 

0 

-20 

-10 

1500 

11 

-30 

6.858 

Std. 

10 

0 

-  5 

1000 

12 

-15 

6.858 

+  10% 

20 

0 

0 

500 

13 

0 

6.858 

-10% 

-10 

0 

5 

500 

14 

15 

6.858 

Std. 

-20 

0 

10 

50 

15 

All  QE's  print  every  two  seconds  plus  impact  at  ground. 

* 

For  all  cases  downwash  is  directed  along  the  local  gravity  vector  regardless 
of  the  aircraft  attitude.  Downwash  is  treated  as  a  step  function  with  the 
listed  velocity  of  downwash  as  amplitude  and  acts  over  a  fixed  distance  of 
6.279  meters. 
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RESIDUALS  IN  BURNOUT  ANGLE 


RESIDUALS  IN  BURNOUT  ANGLE 
(RADIANS) 


RESIDUALS  IN  BURNOUT  ANGLE 


RESIDUALS  IN  BURNOUT  ANGLE 


RESIDUALS  IN  BURNOUT  ANGLE 


tfl 


RANGE  (METERS) 

Figure  5.  Correlation  of  MK  40  trajectory  data  with  single  t2 
term  for  gravity  drop,  wind  not  compensated  for,  and 
coefficient  as  chosen  by  BHT.  (Test  data  superimposed.) 
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RESIDUALS  IN  BURNOUT  ANGLE 
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MS  Amy  Lopi sties  evaluation  Apencv 

nod  Hevclonnent  Command 

ATTM :  OALO-LF.I ,  Mr.  Merritt 

A"1’'1 ;  m»nr:o-pPA-SA 

Fort  Monmouth,  N.T  07701 

Mew  Cumberland,  RA  17070 

1 

director 

Commander 

US  Arnv  TRAOOC  Systems 

"S  Armv  electronics  Research 

Analysis  Activity 

and  Oevelopnent  Command 

ATTN:  ATAA-SL ,  Tech  Lib 

Technical  Sunnort  Actlvitv 

'.'bite  Sands  Missile  Ranp.e 

ATT.':  OrLSn-L 

Fort  "onnouth ,  H.l  07  701 

NM  RR002 

Commander 

O 

Commander 

MS  Amy  Aviation  Center 

MR  Armv  Missile  Command 

and  School 

AT'^tj  0RS'<I-R 

A?mN;  ATZO-n-MS,  Mr.  Warner 

pn^MT-Pf 

ATZO-n-MT,  Mr.  Carbarini 

npSMT-l’np 

Redstone  , Arsenal,  AT,  35898 

’’’ort  Rucker,  AL  Ifi3ft2 

1 'wmmnpdcr 

1 

Commander 

Naval  Orriance  station 

MR  Army  "issile  Command 

ATT'l:  Code  A21,  Mr.  Rerpbnuer 

ATT' :  nP?"T-l!7,  Mr.  Re  rprum 

Ms.  Hrock 

Indian  "ead,  >,T>  20ft/(0 

Redstone  Arsenal,  AT,  35898 

1 

AFF,LM,  The  Rand  Corpnratior. 

ATTM:  Llbi'arv-I) 

Commander 

1700  Main  Street 

HR  Arnv  ''issile  Command 

ATT»?j  nnRMT-YPL 

Santa  Monica,  CA  0040ft 

Redstone  Arsenal,  AL  35898 

n 

Hell  Helicopter  ^extron 

HMT  -  Proiect  Fnpineer  AM-1S 

Commander 

ATTM:  Mr.  Milton  Lamm 

MR  Arnv  Tank  Automotive  Mach 

Mr.  Ceorpe  Cohill 

and  Oevelopnent  Connand 

P.0.  Hox  412 

ATT':  pPPTA-ML 
"nrren,  VI  Asoop 

Ft.  Worth,  t::  7 ft  1 0 1 

DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

3  Hughes  Helicopters 

ATTN:  Mr.  Donald  Price 
Mr.  Dean  Zes 
Mr.  R.  Watson 
Bldg.  305  M/S-T73A 
Centinella  Teale  St. 
Culver  City ,  CA  90230 

2  Teledyne  Corp. 

ATTN:  Mr,  Roger  Gilreath 
Mr,  Tor  Kinsinger 
19601  Nordhoff  St. 
Northridge,  CA  9132A 


ABERDEEN  PROVING  GROUND 

Dir,  USAMSAA 

ATTN :  DRXSY-D 

DRXSY-MP,  n.  Cohen 
DRXSY-AAG,  Mr.  Lauzau 
Mr.  Nicholson 
DRXSY-J ,  Mr.  O'Bryon 
Cdr ,  USAfECOM 

ATTN:  LRSTE-TO-F 

DRSTE-CT-A,  MAJ  C.  Borns 
Mr.  J.  Farrell 
Mr.  J.  Lee 

Dir,  USACSL,  Bldg  E3516,  EA 
ATTN :  DRDAR-CLB-PA 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.)  _ _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


